Resonances in Extreme Mass-Ratio Inspirals: 
Asymptotic and Hyperasymptotic Analysis 
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An expected source of gravitational waves for future detectors in space are the inspirals of small 
compact objects into much more massive black holes. These sources have the potential to provide a 
wealth of information about astronomy and fundamental physics. On short timescales the orbit of 
the small object is approximately geodesic. Generic geodesies for a Kerr black hole spacetime have 
a complete set of integrals and can be characterized by three frequencies of the motion. Over the 
course of an inspiral, a typical system will pass through resonances where two of these frequencies 
become commensurate. The effect of the resonance will be to alter significantly the rate of inspiral for 
the duration of the resonance. Understanding the impact of these resonances on gravitational wave 
phasing is important to detect and exploit these signals for astrophysics and fundamental physics. 
Two differential equations that might describe the passage of an inspiral through such a resonance 
are investigated. These diff'er depending on whether it is the phase or the frequency components of a 
Fourier expansion of the motion that are taken to be continuous through the resonance. Asymptotic 
and hyperasymptotic analysis are used to find the late-time analytic behavior of the solution for a 
system that has passed through a resonance. Linearly growing (weak resonances) or linearly decaying 
(strong resonances) solutions are found depending on the strength of the resonance. In the weak- 
resonance case, frequency resonances leave an imprint (a resonant memory) on the gravitational 
frequency evolution. The transition between weak and strong resonances is characterized by a 
square-root singularity, and as one approaches this transition from above, the solutions to the 
frequency resonance equation bunch up into families exponentially fast. 
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PACS numbers: 04.30.-w,04.50.Kd,04.25.-g,04.25.Nx 



I. INTRODUCTION 



Divergent series are the invention of the devil, and 
it is shameful to base on them any demonstration what- 
soever.^^ Niels Hcndrik Abel's 1828 statement [l| sug- 
gests that asymptotic analysis, which commonly leads 
to divergent series, should not be applied to problems 
of physical interest. Asymptotics, however, has become 
an invaluable tool for physicists seeking approximate an- 
alytic solutions. Multiple-scale analysis, which includes 
boundary-layer theory and WKB theory, allows us to un- 
derstand diverse problems, such as semiclassical quantum 
theory, airplane wing design and turbines Q • In the con- 
text of general relativity, asymptotic (post-Newtonian) 
series [1^31 constitute the basis of the filters used in cur- 
rent gravitational wave detectors to extract signals from 
the noise. 

Resonances are a common occurrence in physical phe- 
nomena. In a traditional oscillatory system, a resonance 
is a point in frequency space where the system stores and 
transfers energy between kinetic and potential modes, al- 
lowing a small driving force to generate large amplitude 
oscillations. Resonant phenomena can occur in many vi- 
brational or wave-like systems and include electromag- 
netic resonances, nuclear magnetic resonance, electron 
spin resonance, and so on. In the realm of general rela- 
tivity, black holes can sometimes be treated as resonators. 



as they relax after being perturbed (see, for example, [H). 

In the context of general-relativistic orbital mechan- 
ics, resonance has recently been adopted to represent a 
slightly different phenomenon, namely, the enhancement 
of gravitational-wave energy dissipation due to the lack 
of cancellation of oscillatory modes that for generic in- 
spirals average out d, 0| • This is particularly relevant for 
extreme mass-ratio inspirals (EMRIs), in which a small 
compact object, such as a stellar- mass black hole or a 
neutron star, orbits around a supermassive black hole [1]. 
In such a two-body system the smaller object slowly spi- 
rals inwards due to gravitational- wave energy-momentum 
losses, on a radiation-reaction timescale much longer 
than the orbital one. This inspiral is usually modeled 
by computing an orbit-average of the gravitational- wave 
energy flux. This procedure discards terms proportional 
to odd-powers of sines or cosines of the sum of the orbital 
phases [9Hl2[ . For orbital configurations or points in fre- 
quency space at which the sum of the orbital frequencies 
vanishes, the orbit-averaged energy flux is not equal to 
the limit of the orbit-averaged fluxes for nearby, nonres- 
onant orbits. This is because harmonics of the frequency 
that vanish on resonance contribute to the secular com- 
ponent of the change in the orbital elements on resonance 
0, 0], but average to zero for off-resonance orbits. Un- 
like traditional oscillators, however, there is no external 
driving force in the EMRI case; the emitted gravitational 
waves drive the inspiral themselves and the resonance is 
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caused by the orbital frequencies becoming commensu- 
rate, that is, some Hnear combination of the three fre- 
quencies with integer coefficients vanishes at resonance. 

A secularly growing radiation-reaction force can leave 
strong imprints on the orbital motion, even if this secu- 
lar growth is active for a very short time. These imprints 
can then propagate into the gravitational waves emitted 
and could have important consequences for gravitational- 
wave detection. Unlike conventional telescopes operating 
in the electromagnetic spectrum, current gravitational- 
wave detectors will not observe signals above the average 
noise. Instead, signals arc expected to be buried deep 
in the noise, and will be extracted using filters based on 
the expected signals. Although EMRIs are not expected 
to be detected with current ground-based gravitational 
wave detectors, they are a key target for future space- 
based detectors for which accurate EMRI filters will be 
needed. It is therefore important to understand how 
EMRI resonances can affect the emitted gravitational 
waves. 



Resonances in Extreme Mass-Ratio Inspirals 

The extrcmc-mass-ratio (typically 10~® - 10"'"') en- 
sures that over short timescales the orbit of the smaller 
object in an EMRI system is approximately geodesic. It 
is therefore appropriate to use an "osculating-element" 
formulation in which the EMRI is identified by a sequence 
of geodesies 0, [HI . Geodesies in a Kerr background are 
uniquely characterized by three constants of the motion, 
energy E, z-componcnt of angular momentum L^, and 
Carter constant Q, and four initial phases that specify 
the coordinates of the object at a particular time. The 
position and velocity of an object uniquely identifies a 
geodesic. Since the evolution of the orbit is governed by 
a second-order differential equation, the values of these 
seven geodesic parameters at each point on the inspiral 
provide an alternative parametrization of the inspiral. 
The time-evolution of the geodesic constants of the mo- 
tion is 

^ = .^f(,,j) + o(a, (1) 

where = {E,Lz,Q) is a vector of these constants, 
while e is the mass ratio, q is an angle phase variable, 
and J^^^ is the "self- force" . The rate of change of 
can then be used to construct the rate of change of the 
orbital frequencies in a similar form. 

For any given geodesic, the self-force can be expanded 
in a Fourier series in terms of the fundamental orbital 
frequencies. These frequencies can be mapped to the 
geodesic constants of the motion. A geodesic resonance 
occurs when the ratio of the frequencies of the radial 
and vertical motion is a rational number. The third fre- 
quency, that of meridional motion, is not relevant for 
resonances due to the axisymmctry for the background 
Kerr spacctimc. Henceforth, we only consider the depen- 



dence of the self-force on the two frequencies that can 
lead to a resonance. 

Let us then expand the self-force in a two-frequency 
Fourier series, where one of the frequencies (w) ap- 
proaches zero while the other (/) remains finite: 

^ e ^ G£„ cos [{£uj + nf) t] + Hen sin [{£lu + nf) t] , 
e,n 

_ (2) 

where Ggn and Hin are time-independent Fourier coef- 
ficients that depend on the orbital parameters. Clearly, 
when £uj + nf ^ 0, the cosine and sine terms average out 
for sufficiently long integration times. However, at reso- 
nance, where £ui = —nf (which in this case we take to 
be ?i = w = 0), the cosine function goes to unity, leaving 
a sum of secular (zero-frequency) Fourier coefficients. 

Let us now further assume that the Fourier coefficients 
(Gin, Hin) vary smoothly as the resonance is approached, 
such that they can be expanded as their on-resonance 
values plus a correction of 0(e^/^). Such corrections can 
always be made small by choosing a sufficiently small e, 
independent of the magnitude of w. Similarly, corrections 
from other terms of 0{e) on the right side of ([1]) can be 
ignored. 

We are now left with a number of rapidly oscillating 
terms (those with £ljj + nf ^ 0) and also terms that 
slowly oscillate away from resonance and then vanish at 
resonance. The rapidly oscillating terms are less impor- 
tant because they average to zero on a short timescale. 
Changing variables to y = i^/\/e and x = \ftt, we find 
that Eq. ^ becomes 

d-V ^ ( r, nfx\ ■ ( n nfx\ 

— = 2^ Gin cos \^£xy + j +Hin sm \^£xy + -j=- j . 

(3) 

Expanding the sum for the first few (£, n) modes and 
rescaling hy y — y/VGoo, x = V Gqox, we then have 

= 1 + fc cos {x y) , (4) 

where we have dropped the tildes, prime denotes differ- 
entiation with respect to cc, and, in principle, the param- 
eter k = Goi/Goo is known. In most scenarios fc ^ 1, 
but there could be orbits for which k = 0(1). 

In deriving (j4|), we have made several approximations: 
(a) we have ignored the rapidly oscillating e~^/^ terms, 
and thus considered only the n = modes; (b) we have 
ignored a phase constant induced by the i?io term; (c) we 
have considered only the £ = 1 mode because these are 
the dominant ones and are on resonance for the longest 
time. Assumption (a) is justified, given that the rapidly 
oscillating components tend to average out over a suf- 
ficiently long integration time. The relaxation of as- 
sumptions (b) and (c) will be addressed more carefully 
in Sec. El 

In addition to the approximations described above, ^ 
also makes the critical assumption that the self-force can 
be expanded as a Fourier series in the frequency with ar- 
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gument {£u! + nf)t and that the coefficients of this expan- 
sion are continuous at resonance. An alternative way to 
write the same equation off resonance would be as an ex- 
pansion in the phase with argument £(j}+mp, where {(j), '0) 
are angle variables. For a geodesic, the time derivatives 
of {(p,ip) are the frequencies (</)', "0') = (w,/), but if one 
regards these phase angles as fundamental and assumes 
that the coefficients of that expansion are continuous at 
resonance, one ends up with a slightly different differen- 
tial equation: 



whose asymptotic expansion in the limit x — >■ +oo is 



1 + fc cos ( 



which admits the first integral 



(l) + ksm(j) + (j)'{0) 



(5) 



(6) 



where (f>'{0) is an integration constant. 

Equations @ and ([5]) give two alternative descriptions 
of an EMRI resonance, but they are not equivalent. To 
make this clear, we rewrite ([5]) in terms of y: 



y' = 1 + k cos ( / ydx 



(7) 



This equation is equivalent to Q only in the limit 
xy' <C 1. In this paper we seek solutions to these two 
differential equations in the limit x — > +oo. The equa- 
tions are deceptively simple (they are just ordinary dif- 
ferential equations) but due to the nonlinearity, finding 
exact solutions is impossible. 

This paper describes the solution to both the 
"frequency-resonance" differential equation Q and the 
"phase-resonance" differential equation ^ and is organ- 
ised as follows: The leading-order behavior of the solu- 
tion at late times (x — > oo) for both frequency and phase 
resonances are calculated in Sec. |lll The higher-order 
behavior in k of these solutions is given in Sec. IIIII Sec- 
tion IIVI describes the qualitative change in behavior as 
k transitions from fc > 1 to < 1. Section IVl discusses 
generalisations of the resonance equations and explains 
how the solutions are modified. Section IVII gives some 
conclusions and describes possible future work. 



II. LEADING-ORDER ASYMPTOTIC 
BEHAVIOR FOR LARGE x 

Equations (|4|) and ([7]) describe simple models of nonlin- 
ear resonant behavior and similar versions of these equa- 
tions have been studied before. In fact, equations of the 
form y' ~ f{ax + fHy) or y' = f{y/x) have solutions in 
quadrature. However, equations of the form y' = f{xy) 
cannot be solved exactly. Instead, one relies on asymp- 
totic techniques to understand their behavior. 

The prototypical equation to study with these tools is 



i 



a 

X 



{x +oo), 



(9) 



where a = (n + l/2)7r and n is an integer. For a slowly 
varying solution, as x — > +oo, <C 1, which implies that 
cosxy <^ 1 and thus xy = a ^ (n-|-l/2)7r. In fact, one can 
show that corrections to this asymptotic solution scale 
with powers of (l/x)™ for m > 1 Similar techniques 
can be used to show that the solution to equation y' = 
tan 2xy also behaves as y ~ (2n -|- l)Tr/x as x — )- +00 

MM- 



A. Frequency Resonances 

Let us first consider the case fc > 1 of ([4]). For slowly 
varying solutions y' <C 1, 1 + cosxy <C 1 and then 

arccos(— 1/fc) , , , 

y^ ^ — - as x^oo,fc>l. (10) 

X 

When fc < 1, the above solution does not exist and more 
subtle asymptotic techniques must be used. 

Next, we consider the fc < 1 case. The form of (|4]) 
suggests that y ^ y^ = ax as x ^ +00. Let us then try 
this ansatz, which when inserted in Q gives 



a = 1 + cos(ax^) 



(11) 



Clearly, the j/c ansatz is not a proper solution. We can 
understand this by averaging the cosine term over all x 



a = 1 + (cos(ax^)) , 



(12) 



where the angle brackets stand for averaging. The Fresnel 
cosine function is defined via the integral 



C{x) 



2a 



1/2 



cos(ax'^)c?x' 



(13) 



y = cos xy , 



(8) 



and as X — > +00, C(x) ^ l/2 + sin(7rx^/2)/(7rx). At very 
large x then, 

, 2\ 1 l/27r\^^^ sin(ax^) ^, 

cos(ax^)dx - — H ^ -. (14) 

4 V a / 2ax ^ ' 

Note, however, that the second term still depends on x, 
so our ansatz yc is still not a valid solution. 

These considerations motivate the improved ansatz 
y ^ yi as X — > -I-cxd, where 

yi = (1 + c)x + b+— sin[(l + c)x^ + bx] (15) 

X 

with constants ai, 6, c to be determined. Inserting this 
ansatz into ([4]), we obtain 

c + 2(1 + c)ai cos[(l -I- c)x^ ^bx\+0 (x^^) (16) 
^ fccos[(l + c)x^ 4- bx] cos{ai sin[(l + c)x^ -f fox)]} 
— fcsin[(l + c)x^ -I- fox] sinjoi sin[(l + c)x^ -I- fox]} , 
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where we have expanded the cosine function with stan- 
dard trigonometric identities. Comparing terms of the 
left and right sides of this equation, we see that ai = 
0{k). 

Let us now require that /s <C 1. Since the sine and 
cosine functions are bounded by unity and since ai = 
0{k) <g; 1, we know that ai sin[(l + c)x^ + bx] 1 for 
all X, and we can therefore expand the cosine and sine 
functions on the righthand side of (|16p . Performing the 
expansion, we find that p6|) becomes 

c + 2(1 + c)ai cos[(l + c)x'^ + bx] ^ k cos[(l + c)x'^ + bx] 
-aiksin^[{l + c)x'^ +bx] (17) 

to 0{x~^, k^) from which we infer that ai ^ fc/[2(l + c)]. 
We are then left with 



aik 



{l ~ cos[2{l + c)x'^ + bx]} 



(18) 



which implies that c ^ — fcai/2, and thus, c ^ — 
Substituting this back into ai, we find that ai ~ k/2. 
The second term of (|18p is not included in c because it 
must be canceled by terms of O(fc^) in y, which we have 
neglected here. Our solution to (4) then becomes 



yi 











1 x-\ sm 






' 2x 





, (19) 



with remainders of 0(x~^ , k^), and where b remains un- 
determined and depends on the initial conditions. We 
have solved ^ numerically in the range k G (0, 0.5) 
and X € (0, lO'^) and verified that indeed is a good 
approximation to the numerical solution, as we show in 
Sec. [ml 

The frequency evolution described by (|19p is particu- 
larly interesting. At late times, the behavior of the fre- 
quency is dominated by the term linearly proportional 
to X, with all others becoming subdominant. The slope 
of the frequency, however, is dependent on k. That is, 
as the physical system goes through a resonance, it ac- 
quires a slope correction that depends on the properties 
of the resonance (that is, on k), a resonant memory of 
sorts. If present in EMRIs, this resonant memory could 
have a large impact on the gravitational wave phase as 
the system traverses a resonance. 



1. Matched Asymptotic Expansion 

The constant b is fixed by the initial condition imposed 
at a; = 0; which requires that a solution be valid in the 
a; ^ 1 limit. Recall that the solution found in p9)) is 
valid in the a; ^ 1 limit and that it diverges as x — 0. 
Let us now look for a solution valid for kx ^ I, with 
/c ^ 1, by using the ansatz 



n=0 



fc"?/„(a;) , as /ca; < 1, fc < 1 , (20) 



where yn{x) are undetermined functions independent of 
k. We could have expanded y in fc(fca;)" instead of in /c", 
but this would lead to more complicated differential equa- 
tions, although the solutions would be the same. Here, 
we choose initial conditions ?;ri(0) = for all n, but the 
extension to more general initial conditions is trivial. 
The zeroth-ordcr solution (n ~ 0) satisfies 



The first-order solution (n 
y'i{x) = cos(a;^) ^ 



yo{x) = 
1) satisfies 



yi{x) 



(21) 




(22) 



To next order in fc, we substitute the solutions found thus 
far into the differential equation 



ky'2 



kxyi{x) + k'^xy2{x)) . (23) 



We can expand the cosine using the assumption kx <^ 1 
to find that 



y'2ix) = -xsm{x^)yi , 
which then leads to the solution 



(24) 



(25) 

Similarly, the equation satisfied by the third-order term 
in the expansion is 

^2 _ 

y'six) = — - cos{x'^)yf{x) - xsm{x'^)y2{x) , (26) 

but this cannot be explicitly integrated. Putting together 
all the pieces found so far, we get 




+ k^ 



C0S(X2)C(4= 



-C 



2x 

7^ 



(27) 

Let us now asymptotically match this solution to the 
one found in ([T9)) . For such a procedure to be valid, a 
buffer zone must exist where both solutions are simulta- 
neously valid. Since (|19|) was found by assuming a; 3> 1, 
while (P7)) assumes that fca; ^ 1, this implies that a buffer 
zone does exist with extension 1 ^ x <C k~^. Asymp- 
totic matching requires that we asymptotically expand 
([T9| in fca; ^ 1 and (|27l) in a; ^ 1, and then set these 
two expansions equal order by order. To leading order, 
we find that 



27r, 



k x + b , 



(28) 



with remainders of relative 0{l/x,kx). This immedi- 
ately leads to 



k. 



(29) 
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B. Phase Resonances 

Let us first consider the solution to the phase resonance 
equation, ([5]), in the k > 1 case. A constant solution as 
X ^ 1 exists provided </)" = 0' = 0, which is satisfied if 

(f> ^ arccos(— 1/fc) , 

(j)'{0) - arccos(l/fc) - y/k^ - 1 - tt . (30) 

The first condition is the same one that xy had to satisfy 
in the frequency resonance case, but the second condition 
now imposes a constraint on the initial conditions. If 
one chooses the initial conditions ^(0) = (/''(O) = 0, the 
second constraint leads to 



\/k'^ - 1 = 2mr - arccos(-l/fc) for n e Z+ (31) 

where arccos(x) is the principal value of the inverse co- 
sine, taking values in the range [0,7r]. Expanding for 
fc 3> 1 we find the approximate solution 



Snn^ + 16n^n^ - 8] . (32) 



TT 1 

kth^ --{l-mT) + - [tt^ 

Evaluating this expression for the first few values of n = 
1,2,.. ., we find that fcth 4.60378, 10.9499, . . ., which 
arc to be compared with the exact numerical solutions to 
(|3T|l . which are fcth = 4.60334, 10.9499, .... We sec that 
the error in the above asymptotic expansion goes roughly 
as 1/fc^. 

Let us now consider the solution to ([5]) in the fc < 1 
case. As before, we concentrate on perturbative solutions 
in fc <C 1. The zeroth-order solution is found by setting 
fc = in ([5|): cj) ^ (f>o = x- /2. The first-order solution in 
fc can be found by postulating that 



00 (x) -I- 01 (x) 
Inserting this into ([5]) we have 

'1 o , 



~ k cos [ —X 



(33) 



(34) 



We assume that (f>i is subdominant relative to (f>Q , and so 
we approximate the argument of the cosine as x'^ /2. We 
can then solve exactly for 0i to find that 

„2 \ 1 



(35) 



where C (x) is the Fresncl cosine function defined in (|13p . 

We can now compare this solution to the one obtained 
for frequency resonances. Differentiating and taking the 
X 3> 1 limit, we find that 



0rfc 



Oix- 



(36) 



Notice that as a; oo this agrees with in func- 
tional form, but not in slope. That is, as the system goes 
through a phase resonance, the slope is not corrected by 
fc. We will find in Sec. IIIII that this remains true as one 
calculates the solution to higher order in fc. Therefore, 
although frequency resonances seem to induce a memory, 
phase resonances do not. 



III. HIGHER-ORDER ASYMPTOTIC 
BEHAVIOR FOR x > 1 AND k < 1 

In this section we consider higher-order solutions in the 
X ^ 1 limit for both phase and frequency resonances, and 
then compare these to numerical solutions. 



A. Frequency Resonances 

In order to obtain a higher-order solution, we must 
construct an ansatz that eliminates the a;-dependent part 
of the right side of ((T8| . We thus pose the ansatz y ^ 
Vasy = yi +112, where yi is given in (dH]), while y2 is 



2/2 



02 



sm{2[{l + c)x^ + bx]} . 



(37) 



Inserting this ansatz into (jH), wc find that 



c + 2{l + c)aiCi + 4(1 + c)a2C2 + Oix-^) 

fccos[(l + c).T^ + bx] cos(aiS'i + 026*2) 

fcsin[(l + c)x'^ + bx] sin(ai5i + a2S2) , (38) 



where we use the notation 



5,1 = sin{n[(l + c)x'^ + 6a;]} , 
C„ EE cos{ri,[(l + c)x^ + bx]} . 



(39) 



As before, we note that Sn and C„ are bounded by unity, 
and since oi = 0{k) we expect that 02 = C'(fc^) or 
smaller. This suggests that we can expand the cosine 
on the right side of the above equation as in Sec. Ill Al to 
obtain 

c+(l+c) (2aiC7i + 402^2) - (1 + C2)+kCi , (40) 

to 0{x^^ , k^). We see then that our previous solution 
still holds: c — —ka\/2 and 2(1 + c)ai = fc, implying 
that c ~ — fc^/4. Wc also see that 4(1 + c)a2 = aifc/2, 
which implies that 02 ~ fc^/[16(l + c)^] or simply that 
0-2 ^ fc^/16 when expanding in fc <C 1. The second-order 
solution therefore becomes 



y2 



16a; 



■ sm 



-4 



26a; 



(41) 



We can obtain the next-order solution by constructing 
the ansatz y j/asy =2/1+2/2+2/3, where ys = ^83 and 
where we assume that 03 = 0{k^). Inserting this into 
O, we get 

c + 2(1 + c)aiCi + 402(1 + c)C2 + 603(1 + c)C3 
fcCi cos(ai5'i + 02^2 + a3'S'3) 
— fcS*! sin(aiS'i + 025*2 + 035*3) . (42) 
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Expanding the above equations in ai <C 1 and 03 ^ 1, 
we find that 

c + 2(1 + c)aiCi + 4a2(l + c)C2 + 603(1 + c)C3 
kai f ka2 kal\ fcai 

-^—+\k^- ^ c^ + —c. 



(^ + ^)C3. 



(43) 



Matching cosine coefficients, this leads to the following 
system of equations: 

kai 

_ 02 
2 

Oi 

2 ' 



2ai(l + c) 
402(1 + c) 
603(1 + 0) 



k I 1 
k'- 



8 



(44) 



which we can solve as an expansion in k to find that 
k fc2 
2 ' "^^^ 16 ' 



96 ' 



k^ 
4 64 



(45) 



Therefore, our solution to third order becomes 



Vaay 



k"^ '\ 
Z lz,4 

4 64 



k 1 

2x 
k^l 
L6' 

fc3 1 . 

sm 

96 a; 



sm 

2 X 



16 X 



2 1 



k 3 4 \ 2 

-Trrk x 

4 64 

4 64^ 



fr2 

1 - - - 

4 64 



6a; 



26a' 



36a; 



(46) 



Notice that this higher-order solution retains the reso- 
nant memory computed in the previous section (that is, 
the k^ correction to the linear-in-a; term, which is domi- 
nant at late times). 

This procedure can be generalized to arbitrary high 
order in A; ^ 1. To leading order in l/x, we make the 
ansatz y ~ 2/(1), where 



y(i) = (l + c)a; + 



n—1 ^ ^ 



(47) 

where (ai"'^'', 61^'', 6, c) are constant coefficients that de- 
pend on k. We can expand {an \ bn '^) as an expansion in 

fc, that is, = X^m^™"^™' ^^'^ solve for these coeffi- 
cients by equating the coefficients of the S'„'s and C„'s 
at different orders in k. Because the derivative of (l47l) is 



2n(l + c)fc 



(aWC„-6«5„)+o(l) 
(48) 



we can evaluate (|4]) to obtain 



2/(1) = I + k Ci cos 



fc^aii)5„ + 6Wc„ 



k Si sin 



(49) 



The generic sine and cosine Taylor expansion formula 
allows us to rewrite the above equation as 



ii,2e 



i-iyk- 



i=0 

(_l)f+ip^+i 



^a«5„ + 6ii)C„ 

00 



2^-1-1 
(50) 



At this point, no further progress can be achieved because 
one needs to evaluate the (2£-|- l)st and the (2£)th power 
of an infinite series, which is not easy to do in closed form. 
This is why it is more convenient to expand the first few 
terms in the series, as done earlier in this section. 

We can generalize the previous procedure to higher 
order in x. Postulate the ansatz y y^x) + J/(2)j where 



yi2)-^Y.^^n^Sn+bi?C^. 



(51) 



The derivative of y is then simply 
2n{l + c)k 



y' = l + c + Y, 

n=l 

2n{l + c)biP + aiP 



00 



(aWC„-6W5„) (52) 



2n(l + c)a^'^ - bi^^ 



while the right side of (|1]) implies that 



y' = 1 + k Ci cos 



k Si sin 
k^ 



X 



-Ci sin 



kJ2{al'^Sn + blpC, 
n=l 

fc^(aW^„ + 6«C, 

00 

fc^(a«5„ + 6WC, 



(53) 



71=1 



n=l 



X 



-Si cos 



00 



n=l 
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One could now expand these equations in fc ^ 1 and 

(2) 

equate coefficients to get equations for tlie a„ 's and 

(2) 

bn s. Following this scheme, one can find the subdoni- 
inant terms in the asymptotic expansion of the solution 
as series in l/x. 



B. Phase Resonances 

Let us now concentrate on higher-order solutions to 
the phase-resonance equation in the fc < 1 case. Wc thus 
postulate the ansatz (j> ^ 0asy , where 

(/-asy = </'0 (x) + (/)! (x) + cj)2 {x) , (54) 

where we recall that 4>q = and we rewrite 0i as 



2 



k 



C 



X 



sm 



(55) 

We have factored out the unbounded-in-x part of 0i and 
the second term is now bounded for all x and tends to 
as a; — )■ 00. This ensures that the term in curly brackets 
is small for sufficiently small fc as x — > oo. We can then 
see that (/)2 must satisfy the differential equation 



k cos ( 



(56) 



where we seek solutions accurate to 0{k^) and, as in 
Sec. IIIB| we neglect the (f)2 term in the source. Inserting 
(/)o and (j)i and using the fact that the bracketed term in 
([55]) is everywhere small to expand the cosine, we find 
that 



k cos 



fc^ sin 



X 

Y 

2 
C 



k^-^x 

-k^x 



(57) 



Integrating this equation twice, imposing the condition 
that 02(0) — 02(0) = 0, and ignoring terms explicitly 
proportional to k^ or higher, we obtain 




The slope of this solution for large x gives the k- 
corrcction to the gradient 



fasy 



1 



k 



Oik') 



(59) 



Notice that although the constant is fc-corrected, the 
linear- in-x term is not, showing again that phase reso- 
nances do not acquire a resonant memory imprint. 

Wc proceed to higher order, and by analogy we write 
down the solution for 02 [x) as the sum of the part on the 
first line of ([55)) that grows linearly with x and a part that 
is bounded for all x and has a convergent integral on the 
range [0,oo]. Schematically, this takes the form 
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2^2 



■+A;$i(x)+fc2$2(a:;), 
(60) 

where fc$i(x) is the term on the second and third lines 
of while k'(^2{x) is the term on the fourth and fifth 
lines. Notice that fc$i(x) does contain terms of O(fc^). 
We now seek the next order solution, 03, that solves the 
equation 



I 1 2 

KCOS -X 

2 



(61) 



where again we have neglected 03 in the source term. 
Inserting the solution known so far and expanding the 
cosine, keeping terms up to O(fc^), we find that 



: k cos ' 

— cos I — x" 
2 V2 

fc^ sin —x' 
\2 



k^ sin I — x^ 
2 



-k 



-k' 



$2 



^kx I $1 



$2 



(62) 



Note that in each term we have eliminated terms in the 
arguments of the cosine and sine that are lower order 
than k^~", where n is the order of the k prefactor to the 
term. By integrating these terms over the range [0,oo], 
we can derive the 0{k^) correction to the asymptotic gra- 
dient. We find the contributions to the O(fc^) term in the 
gradient from each fine of (HH) are 0r/2V2, -0.06202, 
0, and 0F(1 - 2 V2 + V^)/8- 0.07844 respectively. The 
final form for the asymptotic gradient is then 



fasy 



2V2 



0.46484fc^ 



(63) 



C. Comparison to Numerical Results 



The left panel of Fig. [T] shows the difference between 
the numerical solution for y and the asymptotic expan- 
sion in (j46|) to all computed orders in k (solid line) and 
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FIG. 1. Left: Difference between the numerical solution, y, to the frequency resonance equation and the asymptotic expansion 
in (|46p to all computed orders in k (solid line) and to linear order in k (dotted lines) as a function of x. Right: Same difference 
as left-panel but for the solution, </>, to the phase resonance equation. Observe that in both cases the full solution does much 
better than the 0{k) truncation. 



to linear order in k (dotted lines). In all plots we have 
chosen h to be that found in Sec. lII A II As expected, the 
solid lines are much closer to zero than the dotted ones. 
Moreover, notice that the asymptotic solutions found in 
the limit a; » 1 are already quite good (better than 1% 
relative to the numerical solution) below x < 10. 

The right-panel of Fig. [T] shows the same type of dif- 
ference as the left-panel, but for the phase solution. As 
before, observe that the full solution agrees with the nu- 
merical solution much better than its truncated version. 
Unlike the j/asy solution, the (jjusy is globally valid, as we 
did not restrict attention to the a; 3> 1 limit, and this 
can be clearly seen in Fig. [TJ It is also apparent that, as 
expected, in both cases the smaller k is, the better the 
asymptotic solution. 

Although the previous figures establish that the 
asymptotic solutions are indeed accurate representations 
of the numerical ones, they do not compare the y and cj) 
solutions to each other. Figure [5] shows the difference be- 
tween y and d(f)/dx computed numerically in both cases. 
As predicted by the asymptotic solutions, the difference 
is approximately k'^x/2 (the dotted lines) in the a; ^ 1 
limit. 



IV. ASYMPTOTIC TRANSITION AT A: = 1 

As described in the previous sections, both the y and 
(j) solutions experience a transition as fc — > 1. In this 
section we discuss this transition in more detail using 
hyper-asymptotic tools. 
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FIG. 2. Difference between the numerical solution for y and 
the a;-derivative of the numerical solution for tj). For compar- 
ison, we also plot the asymptotic slopes k^x/2. 



from > 1 to fc < 1 from above. To do so, we must find 
the most important term as k approaches 1 to all orders 
in powers of 1/x. Then, we must sum the series and 
identiiy the singularity at fc = 1. 

Assuming that fc > 1, we know that the behavior of 
y{x) in (U) as a; — > CO is described by a series in inverse 
odd powers of x: 



a; x'^ x^ x' 



X — ^ oo. 



(64) 



where 



A. Frequency Resonances 



cos(a) 



(65) 



Let us try to understand the fundamental change in 
the asymptotic behavior of the y solution as k transitions 



As discussed earlier, asfc— )-l+,a— )■ (2ri,-|-l)7r, and thus 
sin a approaches 0. To higher order in 1/a;, one easily 
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finds that 



k sin(a) 

6a sin(a)fc + 



(66) 
(67) 



2fc3[sin(a)]3 ' 
Thus, when sin(a) is small, the most singular part of c is 

2fc3[sin(a)]3" ^ ' 

Similarly, the most singular part of d, e, /, . . . is 
a3 5a^ 



/ 
h 



2fc5[sin(a)]5' 

8fc9[sin(a)]9' 

33a^ 
16fci3[sin(a)]i3 



8F[sin(a)]7' 
16fcii[sin(a)]ii 



The numerical coefficients, 1 
a very simple formula: 

(2n)! 



1 1 

2' 2 



(69) 

. , are given by 



F(n) 



n!(n + 1)12'^ 



{n^O, 1, 2, 3, 



Therefore, if we sum the most singular terms as k 
to all orders in powers of l/x, we find that 



y(x) h 



X fca;3 sin(a) 
From the Taylor expansion 



Fin) 



c^A;2[sin(a)] 



(70) 
1+ 

(71) 



1 



n=0 



one finds that 



F(n)^" 
(72) 



y(2^) - 

X 



a + k sin(a) — k sin(a) W 1 ~ 



2a 



a;2(fc2 _ 1) 

(73) 

where we have used the identity fc2[sin(a)] 2 = fc^ — 1. This 
shows that there is a square-root-branch-cut singularity 
where the asymptotic behavior goes complex as k —> 1+. 

In addition to this branch-cut singularity, one can also 
show that the higher-order terms in l/x bunch up into 
families as a; ^ oo, with each pair of families separated 
by an unstable separatrix curve. As shown earlier, as 
fc — > 1+, the leading-order slope of the solution a — ?> 
{2n + 1)tt and sin a approaches 0. There are, however, 
many solutions for a in (|65l) as k — )■ 1+: the first lies 
just below TT (but above tt/2); the second lies just above 
tt; the third and fourth lie just below and just above Stt, 
and so on. 

Consider two different solutions, yi{x) and 1/2(2;) corre- 
sponding to one of the infinite number of possible values 
of a and define 



Y{x) = yi{x) - V2{x). 



(74) 



Observe that Y{x) satisfies the differential equation 

Y'{x) — k cos [xyi{x)] — k cos [xy2{x)] . (75) 
Using the identity 



COS a — cos p = — 2 sm | — - — | sm 



we can rewrite ([75]) as 

y' {x) ^ — 2fc(sina) sin 



xY{x) 



(76) 



(77) 



for large x. 

Let us now make the assumption that yi [x) and 2/2 {x) 
approach one another gets large, so that Y{x) is 

small when x >> 1. Then, this differential equation 
becomes 



Y'{x) - x(tana)y(a;). 



whose solution is 



Y{x) - C 



(78) 



(79) 



Note that this solution is growing exponentially if tan a 
is positive, and thus the assumption that Y(x) is small 
as X — > C!0 is not valid. This is the unstable (separa- 
trix) case. However, if tana is negative, then we have 
the stable case, and we have shown that the family of 
solutions corresponding to this case all bunch together 
exponentially fast. Note that there is an alternation be- 
tween stable and unstable behavior: Stable behavior oc- 
curs only for the values of a that are just below {2n + l)7r 
for 71 G Z, that is, tt, 37r, Stt, and so on, while unstable 
behavior occurs for the values of a just above (2n -I- l)7r. 



B. Phase Resonances 

The phase solutions also show a fundamental quali- 
tative change of behavior as k transitions from k > 1 
to k < 1. Because the phase differential equation, ([5|), 
admits a first integral of the motion, the analysis is sim- 
pler than in the frequency case and it does not require 
a hyper-asymptotic analysis. This is best understood if 
we consider the first-order form of which we rewrite 
here as 



E^V{(i)\k). 



(80) 



where V{4)]k) = —0 — fcsin^ and E ~ (t>'{^)- The po- 
tential is shown for representative values of k in the left 
panel of Figure |3l Motion can only exist in regions where 
E > V{(f>). The initial conditions (j){0) = 0'(O) = cor- 
respond to E = 0. For fc < 1, the effective potential 
has no turning points and for any initial conditions the 
motion will be unbounded with (f> — >■ +00. When /c > 1, 
the potential does have turning points. The right panel 
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2 4 6 8 10 12 2 2.5 3 3.5 4 

<|, 

FIG. 3. EfTective potential for the phase resonance equation, as defined by Eq. (|80p . The left panel shows the potential for 
four different values of k, while the right panel shows a close up of the region 2 < < 4 when k = 1.1. Each horizontal line 
corresponds to a particular choice of the constant E. Motion can only exist where E > V[4>;k), as indicated by the solid parts 
of the lines shown. 



of Figure [3] shows a close-up of the potential for k = 1.1 
in the vicinity of the turning points. For E > —3.11303, 
the motion is unbounded as in the fc < 1 case. However, 
for —3.11303 > E > —3.17015, the motion intersects the 
potential twice, and it can be oscillatory for suitable ini- 
tial values of (p, or unbounded if 0(0) is sufficiently large. 
For E < —3.17015, the motion is again only unbounded 
for sufficiently large 0(0). 

These specifications for E and 0(0) place restrictions 
on the initial conditions, which are inconsistent with the 
conditions we want to impose. For the conditions 0(0) = 
0'(O) = 0, the motion is unbounded for k < 4.60334, the 
critical value computed in ([5^ . When k « 4.60334, the 
potential intersects the E = axis a second time. The 
motion will asymptotically approach the limiting value 
0CX3 « 4.49341. For k > 4.60334, the motion is oscillatory. 
At the next limiting solution for k, k ^ 10.9499, the 
effective potential has another intersection with the E = 

axis. However, this region is inaccessible to motion with 
these initial conditions and the motion is still oscillatory. 
Starting the motion with 0'(O) = and 5.73224 < 0(0) < 
10.9041 would generate a solution that asymptotically 
approaches the next limiting value 0oo ~ 10.9041. 

At the threshold value k « 4.60334, two solutions 

01 (2;), 02 (a^) with > El ^ E2 < will oscillate with 
different frequencies, and we therefore expect the differ- 
ence 4>i{x) — 02(a;) to be oscillatory. If iJi < < E2, 
one solution will be unbounded and therefore the differ- 
ence will grow like For < Ei ^ E2 > 0, both 
solutions are unbounded, and we expect the difference to 
grow linearly with x. 



V. GENERALIZED RESONANCES 

Section |T] derived certain equations [Eqs. @ and ([5])] 
that are representative of phase and frequency resonances 



in EMRIs, but in doing so we made two important sim- 
plifying assumptions: (b) we ignored the sine term in the 
sum given in Eq. ([3]) and (c) we ignored higher-^ terms 
in this same sum. In this section we relax these two as- 
sumptions and discuss how the solutions are modified. 

Let us first relax assumption (b). If the sine term is 
included, ([4]) and ([5]) can be written as 



4^ = 1 + fccos(xy -f (50), = 1 

ax dx^ 



kcos{4> + (50) , 
(81) 

where Scf) is a constant. We can repeat the analysis of 
the frequency evolution equation with the modification 
introduced above by making the ansatz 

yi = {l + c}x + b+— sin[(l + c)x'^ + bx + <50] (82) 

X 

in place of The analysis proceeds exactly as before, 
but with the arguments of the various cosine and sine 
terms modified via (l+c).T^+6a; —> {l+c)x'^+bx+5cj). The 
asymptotic slope, (l + c), is unchanged as a function of k. 
The solution for b will be modified, however, because the 
expansion for kx 1 described in Sec. Ill AT] is modified. 
In particular 

yo{x) = x, 

1 ^ 

y2{x) 



in which S{-) denotes the Fresnel sine function, defined 
by a similar equation to Eq. (|13p . but with the cosine 
replace by a sine. After asymptotic matching, we then 
find that b \/7r(fc/2) cos(7r/4 + S<j)). 




11 



In the phase-resonance case, the addition of the S4> 
to the equation of motion is equivalent to solving the 
original problem with a modified initial condition: (t>(0) = 
S(j> and (p'{0) = 0. However, the solution to the modified 
equation with the standard initial condition 0(0) = 
can also be found straightforwardly using the method 
described in this paper. In that case, the zeroth-order-in- 
k solution is unchanged, (/<o = 2;^/2, but the first-order- 
in-fc correction, (Ai , is modified to 



-fcsin 



cos{5(f>)C 

2 



sm{S(l>)S 



64>] +k sm{S(f)) 



(84) 



from which we see that the asymptotic correction to the 
gradient is modified to k^/n/2cos{6(j) + 7r/4). Continu- 
ing to the next order in k, we obtain for the asymptotic 
gradient 



dx 



k\ — cos 
2 



0? 



+ (2 - \/2) cos (^2(50 + 



(85) 



for x 1. We note that the leading-order- in-fe correction 
to the asymptotic gradient in the phase equation and the 
leading-order correction to h in the frequency equation 
can be made to vanish for dcj) = it /A. However, even with 
this choice, the next-to-leading-order correction does not 
vanish, and so although the effect of the resonance can 
be suppressed for certain values of 5(j), it cannot be elim- 
inated. 

Let us now relax assumption (c); that is, let us in- 
clude higher i modes in the resonant differential equa- 
tions. Consider first an evolution equation of the form 



dj/ 
dec 



1 -|- kcos{ixy) 



where n is some integer. By writing Y 
we can rewrite this equation as 



dY 
dX 



= l + kcos{XY), 



(86) 

Vly, X = Vex, 



(87) 



reducing it to the same form that we considered before. 
We note in particular that since the solution to this equa- 
tion behaves as y ~ {l + c)X for x ^ 1, this implies that 
y ^ {I + c)x and hence the value of c is unchanged as a 
function of k. The value of b would be modified, however. 
In the phase case, if we modify the equation to 



dx^ 



1 + fc cos 



with n again an integer, the substitution $ — 
\/ex gives 



d2$ 
dJO 



1 + A:cos($). 



b,X^ 



(89) 



Using the solution to this equation that we found earlier, 
we find that (f> w X^/2 + a{k)X + ■ ■ ■ and we deduce 
that 



1 



a(k) 

~7T 



X 



(90) 



Clearly then, higher-^ modes in the sum of ^ are sup- 
pressed by a factor of 1/ VI. 

Let us now consider the case where we have more than 
one term on the righthand side of the evolution equa- 
tion. In the frequency resonance case, we would have an 
equation like 



dy 
dx 



— ~ I + ki cos{xy) + ki cos(f xy). 



(91) 



To solve this equation, we can proceed as before, making 
an ansatz of the form 

j/i EE (1 + c)x + 6 + sin[(l + c)x^ + bx] 



^ sm[e{{l + c)x'^ + bx)]. 



(92) 



Inserting this into the differential equation, we find the 
same source terms involving sin^(-) that we found earlier, 
one from each mode. We also find various cross source 
terms, but these do not contribute to the zero-frequency 
part of the solution because they are products of oscil- 
latory functions with unequal frequencies. We deduce 
then that the solution for c at leading order is the sum 
of the solutions treating each of the modes individually. 
At the next order, cross terms will come in that may be 
important, but these will be subdominant and further 
exploration of these is beyond the scope of this paper. 
In the case of the phase equation we would have 



dV 
da;2 



1 + ki cos(0) -I- ki cos 



(93) 



The linear-in-fc term in the solution of this equation can 
readily be seen to be the sum of the linear-in-A: solutions 
to the equation with only one of the cosine terms on 
the righthand side. Cross terms again come in at higher 
order, but by the preceding argument the size of the cor- 
rections from the ke term are suppressed relative to the 
dominant mode by factors of 1/ VI. 



VI. CONCLUSIONS 

We have studied the behavior of the solution to two dif- 
ferential equations that describe the gravitational-wave 
phase and frequency evolution during an EMRI that ex- 
periences a resonant transition. We have found two gen- 
eral differential equations that might describe this behav- 
ior: one based on the assumption that the coefficients 
of a Fourier expansion of the self-force in frequency are 
continuous at a resonance, and the other based on the 
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assumption that it is the Fourier phase coefficients that 
are continuous at resonance. We have solved both differ- 
ential equations at late times using asymptotic methods. 

Depending on the strength of the resonance (controlled 
by the parameter fc), both resonant equations lead to so- 
lutions with qualitative different behavior: Weak reso- 
nances {k < 1) lead to linear temporal growth of the fre- 
quency and quadratic growth of the orbital phase; strong 
resonances (fc > 1) lead to linear temporal decay of the 
frequency function, leading to a constant-phase offset at 
late times. 

Even though the differential equations for phase and 
frequency resonances might both describe the evolution 
of an EMRI through a resonance, we have found that 
the evolution depends on which differential equation one 
assumes. That is, at late times and in the weak reso- 
nance case, the evolution of the frequency in the phase 
and frequency resonant cases possess different asymptotic 
slopes. The difference in slope depends on the value of 
fc, with the frequency evolution acquiring a fc-dependent 
memory in the frequency resonant case that is absent for 
phase resonances. Further work is required to explore 
which of the two equations is in fact most applicable to 
the EMRI resonance problem. 

We also studied the transition between weak and 
strong resonances. We found that at the transition point 
fc = 1, there is a square- root branch cut in the solution 
to the frequency resonance equation. Close to this point, 
we proved that frequency solutions bunch up into fami- 
lies that decay as 1/a; exponentially fast. In fact, there 
is an alternation between stable (bunching up of solu- 
tions) and unstable behavior, depending on the branch 
of solutions considered. 

Future work should concentrate on exploring which of 



these equations is applicable to EMRI evolutions in prac- 
tice and what the implications arc for the construction 
of waveform template models of EMRI signals. The ex- 
istence of a memory effect in the frequency resonances 
is particularly interesting and would have a profound 
impact on our ability to detect EMRI signals. An ap- 
proximate post-Newtonian prescription for the self-force 
on resonance has been suggested [6(] and would provide 
a suitable framework in which to explore these ques- 
tions further. Whichever of the two equations applies 
to the problem that motivated this work, the results de- 
scribed in this paper provide important insights into the 
behaviour of the solutions to these differential equations 
and predictions for the change in the frequency and phase 
of the evolution as the orbit passes through a resonance. 
These results will be invaluable for constructing approx- 
imate models to describe the evolution of EMRI orbits. 
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